RD-fll££  6M 
UNCLASSIFIED 


DRVTIHE  AND  NIOHTTIHE  NEUTRAL  NIND  AND  TENPERRTURE 
NEASURENENTS  FROH  INCO. .  <U>  AIR  FORCE  INST  OF  TECH 
HRIOHT-PATTERSON  AFB  OH  L  0  BELSNA  DEC  85 
RFIT/'CI/'NR-8S-33T  F/G  4/1 


1/1 


NL 


security  cl  ASSiriCATION  OF  THIS  PACE  flWien  Dal*  Enlerrd) 


REPORT  DOCUMENTATION  PAGE 


t.  R£POR  T  NUMiiER 

AFIT/CI/NR  86-  33T 


2.  GOVT  ACCESSIO 


READ  INSTRUCTIONS 
BEFORE  COMI’LETINC.  FORM 


.  RECIPIENT'S  CATALOG  NUMBER 


4  Title  (tutii  Stihatip) 

Daytime  And  Nighttime  Neutral  Wind  And 

Temperature  Measurements  From  Incoherent 

5.  TYPE  OF  REPORT  A  PERIOD  COVERED 

THESIS/Wj550;A;/0N 

Scatter  Radar  At  300KM  Over  Arecibo 

6.  PERFORMING  O^G.  REPORT  NUMBER 

7.  AUTHORfs; 

Leslie  0.  Belsma 

a.  CONTRACT  OR  GRANT  NUMBERfs; 

9  performiiig  organization  name  and  address 

AFIT  STUDENT  AT: 

University  of  Michigan 

10.  program  ELEMENT.  PROJECT,  TASK 
AREA  •  WORK  UNIT  NUMBERS 

_ 

M.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

AFIT/NR 

12.  report  date 

1985 

WPAFB  OH  45433-6583 

13.  NUMBER  OF  PAGES 

_ 35 _ 

14  monitoring  agency  name  a  AOORESSfi/  dHt«rent  from  CortfroUIttg  Offfco) 

15.  SECURITY  CLASS,  (ol  thia  report) 

UNCLASS 

tSa.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 

16.  01 STRIGUTI ON  STATEMENT  ('o^  Ihia  Report) 

APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 

17.  DISTRIBUTION  STATEMENT  (ol  (he  abalraci  enlarsif  In  Block  70,  U  dllloteni  Item  Kepoil) 

_ ^  ^ _ /  T  f) _ 

10.  SUPPL  EMEN  T  ARY  NOTES  ( 

APPROVED  FOR  PUBLIC  RELEASE:  lAW  AFR  190-1 

Bean  for  Research  and 
Professional  Development 
AFIT/NR,  WPAFB  OH  45433-6583 

19  KEY  WORDS  (Coitfltuie  on  revvrao  aide  if  r^eceaaory  nnd  Identify  by  block  number) 

DTIC 

ont  FlLt  COP^ 


DD  1473  EDITION  OF  I  NOV  65  IS  OBSOLETE 


security  classification  of  this  page  (llJirn  D«l*  Enltttd) 


•J- 

/n-  n' 


kv  '  \?il'  rr*  '  y  ' 

>1^  "*  H  tj  *  .»  .  •V'’  ,  '/i-  -  .  . 

>  .'-V-  .'-V-".''-'.-.  V  -  ‘.'-V  '  *  '  -  ■<  "■ 


ABSTRACT 


~  i 

Areciba's  twenty  acre  epherical  radar  reflector  detects  the 
weak  electromagnetic  energy  reflected  when  pulses  are  transmitted 
into  the  ionosphere.  The  information  contained  in  this  reflected 
energy  reveals  the  electron  temperature,  the  electron/ion  tempe¬ 
rature  ratio,  the  line  of  sight  tlasma  velocity  and  the  electron 
density  among  other  parameters.  This  information,  combined  with 
neutral  densities  obtained  from  the  global  thermospheric  model 
based  on  mass  spectrometer  and  incoherent  scatter  data,  liSIS-63, 
is  used  to  determine  the  neutral  temperature  and  winds  in  the  F 
region.  This  study  pertains  to  the  F2  region  at  a  height  of 
300km  over  Arecibo  (18  N  latitude  with  a  magnetic  dip  angle  of  50 
).  A  few  trends  in  previous  data  are  noted  followed  a  des¬ 
cription  of  the  300km  temperature  and  neutral  wind  for  daytime 
and  nighttime  from  12  days  of  data,  sets  of  three  days  from 
January  1984  and  1985,  September  1984,  and  June,  1984.  The 
experimental  neutral  temperature  is  compared  to  the  MSIB  model 
prediction.  The  observed  temperature  closely  agrees  in  phase, 
but  averages  50K  lower  than  the  model.  '^The  observed  winds  are 
compared  to  NCAR's  Thermospheric  Qeneral  Circulation  Model  (TGCM) 
which  has  incorporated  the  effects  of  tideV  propagating  from  the 
lower  atmosphere.  The  observed  winds  dur^g  equinox  agree  in 
phase  with  the  TGCM,  but  are  a  factor  of  thrde  larger  in  ampli¬ 
tude.  The  ion  drag  force  is  calculated,  usahi  to  resolve  the 
pressure  gradient  force,  and  both  are  discussed  in  relation  to 
the  observed  neutral  wind  patterns.  The  pressure  gradient  force 
is  compared  to  that  predicted  by  the  MSIS  model  and  found  to  be 
of  comparable  magnitude,  but  often  dissimilar  in  phase. 


DAYTIME  AND  NIGHTTIME  NEUTRAL  WIND  AND  TEMPERATURE 
MEASUREMENTS  FROM  INCOHERENT  SCATTER  RADAR 


AT  300KM  OVER  ARECIBO 


raqulrafflsnts  for  Mastor  of  Scianca  dagraa  in 
Atmospheric  Science  at  the  University  of  Michigan 


INTRODUCTION 


Th«  radar  datacts  danaity  f luctuationa  which  for  amall  valuaa 
of  tha  Dabya  langth.  Da  (■69(Ta/N>ffl«  whara  Ta  ia  alactron  tampa- 
ratura  and  m  ia  tha  ion  maaa) ,  occur  aa  Joint  longitudinal  oacil- 
lationa  of  both  alactrona  and  Iona.  Whan  Da  ia  much  laaa  than 
tha  wavalangth  of  tha  ion  oecillationa,  alactroatatic  rathar  than 
praaaura  forcaa  dominata  tha  alactrona  and  thay  follow  tha  ion 
motion  (Walkar,  197G) .  Thua,  acattaring  occura  aa  if  tha  plaama 
conaiatad  of  particlaa  with  tha  alactron  tamparatura  and  ion  maaa 
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(Evana,  1974).  Incoharant  acattar  affacta  only  1x10  of  tha 
vartical  radar  baam  and  tha  acattaring  occura  in  all  diractiona; 
howavar  conatructiva  intarfaranca  occura  for  wavaa  propagating 
along  tha  diraction  of  tha  incidant  radiation  with  wavalangtha 
aqual  to  1/2  that  of  tha  incidant  baam,  anabling  tha  acattaring 
phanomanon  to  ba  obaarvad.  Tha  raflactad  aignal  travala  at  tha 
apaad  c,  thua  maaauring  tha  acho  dalay  datarminaa  tha  haight  at 
which  tha  acattaring  occura.  Charactariatica  of  tha  raflactad 
apactrum  ravaal  tha  varioua  F2  paramatara. 

Tha  alactron  tamparatura,  Ta,  ia  datarminad  from  tha  Dopplar 
broadaning  of  tha  apactrum,  which  ia  laaa  than  would  ba  axpactad 
from  random  tharmal  anargy  of  the  alactrona  bacauaa  of  tha  pra- 
aenca  of  iona.  The  difference  in  Te  and  the  ion  tamparatura,  Ti, 
ia  evident  in  tha  apactrum  ahape,  giving  it  a  double  humped 
appearance.  The  plaama  wave  velocity  ia  proportional  to  tha  aum 
of  tha  ion  and  alactron  tamper aturaa,  which  ahifta  the  apactrum 
to  either  aide  of  the  tranamittad  frequency.  However,  damping 
from  the  interaction  of  tha  wave 'a  electric  field  with  tha 
charged  particlaa  decraaaaa  tha  ahift  and  apraada  tha  apactrum 
toward  the  tranamittar  frequency  (Walker,  1978).  Thia  damping  ia 
proportional  to  tha  wave  velocity  and  alao,  therefore,  the 
electron  temperature.  The  damping  ia  weaker  aa  Te  increaaea  and 
tha  apactrum  haa  a  aharpar  peak  toward  tha  ion  wave  velocity  (ia. 
the  wing  of  tha  apactrum  incraaaaa) |  more  damping  occura  for 


lower  Te  and  the  spectrum  is  spread  more  toward  transmitter 
frequency  (Walker,  197B) .  The  shift  is  proportional  to  the  line 
of  sight  plasma  drift  velocity  as  well  as  the  plasma  frequency. 
Therefore,  assuming  the  ions  and  electrons  move  together,  the 
plasma  drift  can  be  measured  from  the  spectrum  shift.  The  power 
of  the  reflected  signal  is  proportional  to  the  electron  density. 


TEMPERATURE  AND  ELECTRON  DENSITY 


PREVIOUS  OBSERVATIONS 

In  their  study  o-f  exospheric  temperatures  over  Millstone  Hill« 
Selah  end  Evans  (1973)  found  that  the  time  of  Tmax  lags  the  time 
of  maximum  heating,  occuring  about  three  hours  before  thermos¬ 
pheric  sunset  and  that  Tmin  occurs  about  two  hours  before  ther¬ 
mospheric  sunrise.  Tmax  occurred  at  i530LT  in  winter,  1600LT  at 
equinox,  and  1700LT  in  summeri  Tmin  occurad  at  0400  in  winter, 
0200  at  equinox,  and  0030  in  summer.  The  amplitude  of  the  daily 
oscillation  exhibited  a  seasonal  variation  with  maximium  ampli¬ 
tude  occuring  in  summer.  They  found  nightime  fluctuations  on 
some  days  occurring  Just  after  midnight  with  a  50K  amplitude  and 
suggested  that  they  could  be  attributed  to  an  over estimation  of 
temperature  resulting  from  the  neglect  of  frictional  heating  due 
to  neutral  wind.  The  seasonal  variation  of  the  diurnal  oscil¬ 
lation  is  not  represented  in  thermospheric  models  and  is  at¬ 
tributed  to  larger  amounts  of  heat  being  deposited  in  the  summer 
auroral  zone  than  the  winter  (Evans  et  al . ,  197B>. 

In  a  study  of  the  night-time  F2  region  over  Arecibo  based  on 
19  nights  of  data  during  sunspot  minimum  years,  the  400km  ion 
temperatures  were  found  to  be  slightly  higher  than  model  pre¬ 
dictions  and  around  midnight  remained  nearly  constant  for  four 
hours  rather  than  exhibiting  the  sinusoidal  variation  of  the 
Jacchia  model  (Fukato  et  al . ,  1979).  Fukato  suggests  the  "flat¬ 

ness"  could  be  due  to  a  semidiurnal  oscillation.  Thermal  equi¬ 
librium  was  established  less  than  one  hour  after  sunset  in  sum¬ 
mer,  but  three  hours  following  winter  sunset.  Electron  tempera¬ 
tures  increased  rapidly  during  summer  sunrise,  but  began  increa¬ 
sing  two  hours  before  sunrise  in  winter.  The  mean  midnight 
temperature  was  700K  in  summmer  and  770K  in  autumn  and  winter. 

In  studying  night-time  conditions  during  equinox  at  350km  over 


Arscibo,  Harper  -found  a  40K  ion  temperature  increase  occuring  in 
the  three  hours  a-fter  midnight  during  a  time  when  neutral  winds 
are  abating  making  -frictional  heating  an  unlikely  explanation 
(Harper  y  1973) .  He  suggested  that  adiabatic  heating  o-f  neutral 
gas  could  cause  the  temperature  increase. 

Three  nights  o-f  data  over  Arecibo  in  August  1982  showed  the 
average  ion  temperature  between  300-400km  to  be  around  9S0K  -from 
2200-0500AST  (Burnside,  1984).  The  study  also  indicated  an 
average  Ti  o-f  875K  -for  three  nights  in  June  1983,  1025K  -for  three 
January  nights  in  1981  (solar  maximum),  and  850K  in  October  1983. 


A  50km  drop  in  the  F2  layer  height  (represented  by  Zmax)  has 
been  observed  a-fter  midnight  -followed  by  a  rise  at  0400AST  at 
Arecibo  (Behnke  and  Harper,  1973).  They  -found  the  maximum  height 
o-f  the  F2  layer  to  occur  in  the  early  evening.  Burnside  indi¬ 
cates  a  seasonal  variation  in  this  pattern  showing  the  post¬ 
midnight  drop/rise  to  be  gradual  in  June  but  sudden  in  Autumn 
(Burnside,  1981).  Harper  reports  that  the  peak  electron  density, 
N,  increases  during  the  midnight  drop  and  mentions  a  similar  Zmax 
descent  pattern  in  the  a-fternoon  (Harper,  1979),  Night-time 
electron  densities  at  the  F2  peak  over  Arecibo  presented  by 
Fukato  exhibit  a  seasonal  variation  in  the  exponential  decay 
which  begins  in  the  late  a-fternoon,  the  rate  being  the  steepest 
in  the  summer  and  much  smaller  in  the  winter  (Fukato  et  el., 
1979).  An  increase  in  N  a-fter  midnight  appeared  as  a  "small 
hump"  in  the  decay  at  0100  in  winter,  but  only  as  a  decrease  in 
the  decay  rate  -from  2200-0100  in  summer.  He  -found  the  expo¬ 
nential  decay  to  continue  even  after  sunrise  in  summer,  but  to 
stop  two  hours  before  sunrise  in  winter  with  N  becoming  constant 
until  sunrise.  His  calculations  of  Zmax  indicated  the  layer 
ascending  after  sunset  with  the  characteristic  descent  beginning 
at  midnight.  Zmax  increased  before  sunrise,  with  the  ascent 
being  concurrent  with  the  time  that  N  ceased  to  decrease. 
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CURRENT  DATA  REDUCTION  AND  MEASUREMENTS 

Rogar  Burns! da  takas  thasa  maasuramants  by  oparating  tha  radar 
in  two  axparlmantal  modas  probing  tha  ionosphara  frotn  140  to 
SAOkflt.  Most  of  tha  obsarving  tima  is  davotad  to  transmitting  a 
long  pulsa  (50km)  with  an  11km  spacing.  Two  haights  of  tha 
"ovarsamplad"  data  ara  than  avaragad  to  giva  23km  spacing  and  a 
rasolution  of  about  50km.  Tha  antanna  is  at  a  15  dagraa  angla 
from  zanith  and  is  rotatad  360  dagraas  for  samplings  at  12-lB 
azimuthal  positions  which  ara  than  avaragad  for  aach  haight  to 
give  the  electron  and  ion  temparatures.  Ion  velocity  is  calcu¬ 
lated  by  assuming  no  spatial  gradient  over  the  horizontal  sam¬ 
pling  area  and  fitting  a  sine  wave  to  the  line  of  sight  velo¬ 
cities  for  the  azimuthal  angles  at  aach  height  to  resolve  the 
horizontal  components  (Burnside*  1984). 

Tha  ACF  data  provide  maasuramants  of  electron  density  with  a 
50km  rasolution*  but  much  more  accurate  maasuramants  are  avai¬ 
lable  from  tha  Barker  coda  data.  After  every  swing  through  tha 
18  azimuthal  positions  a  short  pulse  is  transmitted  and  tha 
Barker  coding  of  tha  pulse  gives  a  haight  rasolution  of  600m. 
Peak  electron  densities  ara  calibrated  by  an  on  site  ionosonde. 
Tha  Barker  measurements  ara  fitted  to  a  Chapman  profile  using  a 
minimum  haight  between  200  and  250  km  and  thasa  peak  electron 
densities  than  compared  to  the  ACF  values  to  eliminate  unrea¬ 
listic  temporal  increases  or  decreases  in  Ne  arising  when  the 
Chapman  profile  does  not  provide  a  good  fit. 


In  the  currant  study*  as  in  previous  axparimants,  several 
assumptions  are  necessary  to  obtain  the  neutral  temperature  at  a 
given  height.  Ions  provide  the  thermal  coupling  between  elec¬ 
trons  and  neutrals.  The  thermal  structure  is  determined  by 
energy  from  photoionization*  thermal  conduction  and  transport. 
The  electron  temperature  increases  as  photoelectrons  transfer 
energy  to  the  electron  gas.  Through  elastic  collisions  between 
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ions  and  alactrons,  the  ion  temperature  Increases  above  that  of 
the  neutrals.  The  difference  (Ti-Tn)  is  usually  less  than  lOOK, 
though  neglect  of  the  densities  of  N  2  and  0  2  leads  to  a  slight 
underestimate  of  (Ti-Tn)  (Balah  and  Evans*  1973).  It  is  assumed 
that  there  is  only  a  single  ion  species  (atomic  oKygen)  present 
in  the  region  being  studied  and  that  the  ions  are  in  thermal 
energy  balance.  In  using  neutral  density  for  a  single  species 
from  the  thermosperic  model*  atomic  oxygen  is  assumed  to  obey  the 
laws  of  diffusive  equilibrium.  The  ACF  data  directly  yield  Te, 
and  Te/Ti|  Ti  is  calculated  and  the  neutral  density  is  obtained 
for  the  corresponding  height  and  time  from  M5IS  given  the  10.2cm 
flux  and  magnetic  Ap  indices.  The  neutral  temperature  can  be 
calculated  with  these  parameters  using  the  thermal  energy  balance 
equation 

(4.B2  X  10^)  N  (Te  -  Ti )  0.21  C03  (Ti  ♦  Tn)^^^ 

(1)  - — - 

Te'^''^  (Ti  -  Tn) 

where  COl  is  the  atomic  oxygen  density.  Next  it  is  assumed  that 
the  temperatures  vary  linearly  with  height  and  using  a  linear  fit 
for  Te*  Ti *  and  Tn  for  seven  heights  in  January  and  June  of  19B4 
from  2B5-424km  and  for  ten  heights  in  Beptember  19B4  and  January 
19B5  from  262-470km*  the  three  temperatures  for  the  given  alti¬ 
tude  of  300km  are  calculated. 

The  neutral  temperatures  from  this  study  are  shown  in  Figures 
lA-D.  Mean  temperatures  in  the  current  data  set  are  considerably 
less  than  those  calculated  at  Arecibo  during  and  in  the  years 
Just  after  solar  maximum.  Burnside  (19B4)  found  the  highest 
average  night-time  temperatures  in  June;  the  same  is  true  of 
this  study.  The  latest  MSIS  model  incorporates  the  "flatness" 
which  Fukato  demonstrated  and  though  it  matches  the  data  more 
closely  than  a  sinsusoidal  variation  of  earlier  models  would*  it 


does  not  reflect  small  amplitude*  short  time  scale  oscillations 


which  occur  ovornight  at  Aracibo.  Harpar  notad  a  40K  tamparatura 
incraaae  batwaan  midnight  and  3am{  a  aimilar  pattarn  occurs  in 
tha  data,  only  3  hours  aarliar. 

In  all  four  sats  o-f  data  in  tha  currant  study,  ion  tampara- 
turas  raach  paak  valuas  batwaan  1700-1B00A8T.  Tha  dif-faranca 
batwaan  Ti  and  Ta  navar  axcaads  lOOK.  In  Juna  and  Baptambar ,  tha 
ions  and  nautrals  ara  in  tharmal  aquilibrium  from  1900-0500AST, 
whila  in  January  aquilibrium  is  approached  more  gradually  and 
lasts  -from  2100  to  only  0300.  Tha  neutral  temperature  pattarn 
reflected  in  tha  data  agrees  well  with  liSIS  predictions  though 
tha  timing  of  the  peak  temperature  often  varies  by  an  hour  or 
two.  Tha  model  consistently  predicts  temperatures  40  to  BOK 
higher  than  observed  with  the  largest  differences  occurring  in 
the  afternoon.  After  midnight  the  model  indicates  a  very  gradual 
temperature  decrease  to  a  minimum,  but  temperature  maxima  inter¬ 
rupt  this  period  of  cooling  in  the  data  sets. 

The  overnight  equilibrium  temperature  in  January  19B4  averages 
72SK,  then  Ti  begins  to  rise  above  Tn  by  0300AST.  Again,  the 
evening  cooling  trend  is  interrupted  by  a  50K  increase  from 
21-2300,  with  a  leveling  off  or  another  peak  at  2am  before  drop¬ 
ping  to  Tmin  at  3am.  Tn  begins  rising  after  0500,  with  a  sharp 
increase  after  9am  to  a  peak  of  900K  by  noon,  not  cooling  until 
after  1700.  Figure  1C  shows  that  NSIS  neutral  temperatures 
average  about  SOK  warmer  with  BOK  differences  at  the  extrema  and 
no  indication  of  the  accelerated  temperature  increase  in  the  late 
morning. 

The  average  night-time  temperature  a  year  later  (and  closer  to 
solar  minimum)  in  January  19BS  is  650K  with  Tn  reaching  a  peak  of 
B25K  at  1700  but  averaging  only  about  750K  in  the  afternoon.  The 
measured  temperatures  average  about  the  same  as  MS IS  from 
1900-0900,  but  the  predicted  temperatures  are  BOK  higher  in  the 
daytime  with  a  peak  at  1600  when  it  is  observed  closer  to 


IBOOAST.  As  in  January  1984,  tha  ion  temparatura  bagins  to 
axcaad  tha  nautral  Just  a'ftar  3am.  A  well  de'finad  tamper atura 
increase  is  not  evident  after  2100AST  and  a  rapid  drop  occurs 
■from  0130-0300,  followed  by  an  increase  from  0300-0500  which  is 
not  reflected  in  the  MSIS  model  (see  Figure  ID). 

The  average  Tn  throughout  the  night  in  the  September  1984  data 
set  is  650K  and  as  in  June,  Ti  is  equal  to  Tn  from  1930-0530. 
Except  from  September  19-20  when  temperatures  drop  to  minimum  at 
midnight  and  begin  rising  steadily  two  hours  later,  the  tempera¬ 
ture  increase  from  21-2300  is  evident,  though  of  less  amplitude 
than  in  June,  and  a  more  gradual  bulge  also  occurs  around  3am 
(see  Figure  18).  During  the  day  Tn  reaches  a  peak  of  B50K  with 
Ti  averaging  about  30K  warmer.  tISIS  predicts  daytime  neutral 
temperatures  40K  higher  than  observed  indicating  maximum  Tn  at 
1700. 

For  three  June  1984  nights,  the  average  night-time  (22-05AST) 
neutral  temperature  at  300km  is  750K,  with  ions  and  neutrals  in 
equilibrium  from  1930-  0530.  Tn  falls  off  quickly  after  the  late 
afternoon  peak,  but  a  50K  increase  disrupts  the  cooling  trend 
from  2100-2300  and  again  from  0200-0400.  During  the  day,  Ti  is 
about  50K  warmer  than  Tn.  An  unusually  rapid  heating  is  evident 
on  June  27  when  Tn  increases  rapidly  from  05-07am  to  850K  and 
975K  at  1700.  This  sharp  increase  does  not  occur  on  the  previous 
morning.  MSIS  indicates  a  slightly  larger  temporal  gradient  in 
the  morning  and  the  phase  of  temperature  oscillation  matches  the 
data.  However,  MSIS  temperatures  average  about  70K  higher  (see 
Figure  ID) . 

Figures  2A-D  show  the  electron  densities  at  the  F2  peak  and 
the  height  of  the  layer  peak.  In  January,  1984,  the  electron 
density  is  the  highest  in  the  morning  when  the  height  of  the 
layer  is  the  lowest  (between  210  and  240km).  The  peak  densities 
occur  in  the  afternoon  around  6pm  in  the  other  three  data  sets. 


In  J«nu*ry,  19B8  and  Saptambar,  1984,  tha  a-ftarnoon  danaity 
maximum  coincidas  with  an  F2  layar  drop.  In  Juna,  tha  layar 
atays  around  300km  throughout  tha  day.  In  all  four  data  aata, 
tha  F2  paak  raachaa  ia  highaat  altituda  around  midnight  and  than 
undargoaa  a  rapid  100km  drop,  followad  by  a  riaa  or  lavaling  off 
by  3am.  Thia  poat  midnight  drop  and  riaa  ia  conaiatant  with 
pravioua  obaarvationa  (Bahnka  and  Harpar,  1973)  though  a  aignifi- 
cant  aaaaonal  variation  ia  not  avidant.  An  abrupt  riaa  and  fall 
of  tha  F2  paak  haight  occura  aftar  midnight  on  Baptambar  19  whan 
tha  magnatic  Ap  indax  Jumpad  to  36  up  from  3  on  tha  pravioua  day. 


Figure  IB.  Average  Tn  lower  than  in  1984.  Mzucinuin  Tn  occurs  2-3  hours  later  than 
the  peak  in  the  MSIS  curve.  Temperature  increase  at  3am. 


Figure  1C.  A  40K  Tn  increase  2-3  hours  before  midnight  and  again  at  3am.  Maximum 
Tn  occurs  an  hour  later  than  model  prediction. 
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Fig\ire  ^A.  Height  and  density 
of  the  F2  peak  from  16-19  Jan, 
1964.  Hipest  densitiy  occuxs 
just  before  noon. 


Figure  2B.  The  same  for  14-17 
Jan,  1985.  Maximum  density  in 
the  afternoon. 


Figure  2C.  "nie  same  for  24-27 
Sept,  1984.  Maximum  peak  den¬ 
sity  in  the  afternoon  corre¬ 
sponding  to  a  collapse  in  the 
layer  height. 


Figure  2D.  The  same  for  24-27 
Jun,  1984. 


NEUTRAL  MIND 


PREVIOUS  OBSERVATIONS 


Amayanc  and  Vasaaur  datarmlnad  day  and  night  nautral  winds 
using  tha  sama  aKparimantal  mathods  as  tha  currant  study  and 
■found  strongar  -flow  to  tha  polas  in  tha  wintar,  tha  affact  of 
transaquatorial  flow  from  summar  to  wintar  hamisphara  (Amayanc, 
1971>.  Thay  maasurad  ion  drifts  and  tamparatura  gradiants  at 
various  F2  lavals  with  tha  St  Santin-Nancay  radar  and,  with 
nautral  dansitias  from  tha  Jacchia  and  Sloway  modal,  thay  datar — 
minad  tha  diffusion  valocity  thus  infarring  the  nautral  valocity. 
Tha  maridional  winds  wars  prasantad  for  two  days  in  1969  and  ona 
in  1968,  charactaristic  of  low  magnatic  activity.  Thay  found 
polaward  winds  from  7am  to  8pm  in  wintar  raaching  a  maKimum 
valocity  of  lOOm/s.  In  summar,  tha  polaward  flow  did  not  axcaad 
SOm/s  and  ravarsad  to  aquatorward  aarliar  in  tha  day  (around 
3pm).  Tha  sama  phasa  occurrad  during  aquinox  as  in  summar,  but 
tha  maximum  polaward  valocity  was  BOm/s.  Tha  southward  flow  in 
wintar  approached  a  ZOOm/s  peak  abruptly  after  tha  evening 
reversal  in  wintar,  but  tha  aquatorward  flow  increased  more 
gradually  to  a  ZOOm/s  peak  in  summer.  The  maximum  aquatorward 
valocity  during  equinox  was  about  IZOm/s. 

Using  the  same  analysis  technique  as  the  current  study,  Behnke 
and  Harper  (197Z)  determined  maridional  neutral  winds  for  daytime 
on  January  8  and  October  31,  1971  and  for  five  winter  nights  at 
Arecibo  in  1970.  Thay  found  poleward  flow  throughout  both  days 
with  maximum  polaward  flow  in  mid-morning.  The  poleward  flow 
maximum  was  140m/s  in  January  and  the  average  around  60m/s.  The 
strongest  polaward  flow  during  equinox  was  60m/s  with  an  average 
of  35m/B.  Assuming  a  steady  state,  Behnke  and  Harper  determined 
the  October  31  pressure  gradient  as  the  force  balancing  the  ion 
drag  term.  Tha  gradient  was  positive  (poleward)  all  day  with  a 
Zcm/sZ  peak  at  10am  and  a  3cm/s2  peak  at  Zpm.  They  explain  the 


IZ 


morning  flow  in  th*  same  direction  as  the  prassurs  gradient  as  a 
result  o-f  electric  -field  e-ffects  where  a  30m/s  poleward  perpendi¬ 
cular  ion  velocity,  V-perp,  adds  40m/s  to  poleward  -flow  and  20m/s 
southward  V-perp  adds  25m/s  southward  component  to  the  neutral 
wind  in  the  a-fternoon.  Looking  at  the  five  winter  nights,  Behnke 
and  Harper  -found  northward  winds  up  to  BOm/s  in  the  early  evening 
reversing  by  10pm  with  an  equatorward  maximum  of  60m/s  two  hours 
later.  Southward  flow  abated  becoming  zero  or  poleward  again  by 
2am  and  averaging  slighly  poleward  through  the  early  morning. 

Harper  determined  the  meridional  neutral  winds  at  Arecibo  for 
four  nights  in  March  1971  from  Spm  to  Sam  (Harper,  1972).  An 
average  poleward  flow  of  35m/s  occurred  in  early  evening  rever¬ 
sing  at  lOpm  as  in  winter.  The  equatorward  winds  climbed  to 
120m/s  around  11pm  and  dropped  to  zero  three  hours  later  becoming 
briefly  northward  but  oscillating  about  zero  for  the  remainder  of 
the  early  morning.  Harper  determined  that  the  F-layer  achieved 
its  maximum  height  coincident  with  the  maximum  equatorward 
neutral  winds. 

Burnside  (1983)  determined  average  wind  contours  from 
2S0-490km  using  Arecibo  incoherent  scatter  data  for  nighttime 
hours  during  solar  maximum.  Averages  for  data  in  December  19B1, 
February  1982,  and  April  1982  indicate  southward  winds  above 
300km  and  weak  northward  below  (an  average  of  nine  nights  in  June 
1980  and  August  1982,  does  not  reflect  this  pattern).  The  equa¬ 
torward  flow  at  the  higher  levels  is  especially  apparent  during  a 
magnetic  event.  The  plot  for  a  magnetically  disturbed  night  on 
28-29  December  shows  southward  velocities  exceeding  200m/s  above 
370km  but  less  than  25m/s  below  300km  (Burnside  et  al . ,  1983). 


CURRENT  DATA  REDUCTION  AND  MEASUREMENTS 


As  prsviously  msntionsd  ths  lins  of  sight  vslocity  Is  obtsinsd 
■from  ths  spsctrum  shift.  Ths  componsnts  srs  rssolvsd  as 
functions  of  sltituds  assuming  no  horizontal  spatial  gradisnt 
ovsr  ths  arsa  covsrsd  as  ths  antsnna  rotatss  and  no  tsmporal 
changss  in  ths  tims  for  ons  rotation  (about  20  minutss).  Si von 
calculatsd  nsutral  tomporaturost  nautral  dansitiss,  t03  and  CN2], 
ars  obtainsd  from  MSIS  and  ths  diffusion  vslocity  in  ths 
dirsction  of  ths  sarth's  magnstic  fisld  is  calculatsd  using  ths 
0-»-  diffusion  aquation  (Burnsids  at  al .  ,  19B3)  t 

(2)  Vdiff  «  (-Da)sinl  Cl/Hp  l/N(dN/dz) 

+  w/Tr<dTr/dz>  +  l/Tp(dTp/dz)  | 

wharsi  Rsducad  Ti  Tr  ■  <Ti+Tn)/2|  Plasma  Ti  Tp  «  (Ti+Ts>/2| 

w  B  0.37(thsrmal  diffusion  factor) | 

Da  «  <3.02xl0^^)Tn''^/  AC03+19. 9CN^3 (Tn^^^) I 
A  -  1.08  -  O.Mlog  Tn  +  0.0045  (log  Tn>^^^ 

Ths  gradisnt  in  slactron  dansity  at  300km  is  obtainsd  from  ths 

Chapman  fits  to  ths  Barksr  cods  dansity  msasuramants  abova  200km. 
Error  occurs  in  dstsrmining  ths  hsight  and  dansity  of  ths  paak 

whsn  ths  laysr  drops  bslow  200km,  but  calculations  of  ths  vsr — 

tical  dsnsity  gradisnt  at  300km  with  ths  F2  psak  far  bslow  should 
not  bs  significantly  affsctsd. 

Ths  componsnt  of  ths  nsutral  wind  along  ths  magnstic  fisld 
lins  is  takan  to  bs  ths  diffsrsncs  bstwssn  ths  ion  parallsl 
vslocity,  Vpa,  and  ths  diffusion  vslocity.  Tharsfors,  in  ths 
horizontal  msridional  dirsction,  nsglscting  ths  vartical  compo¬ 
nsnt  of  ths  nsutral  wind  (Burnsids, 1984) i 


- »  -  > 


■v^rw-S^.y^j  : 


<3>  Ux  -  <Vp«-Vdi-ff  >/cos  I 

The  calculated  neutral  temperature  rather  than  model  calcu¬ 
lations  of  T-in-finity  are  used  to  get  the  most  accurate  neutral 
densities  from  liSIS  (Burnside,  1984).  However,  realizing  pos¬ 
sible  error  due  to  local  variations  in  density  which  are  not 
reflected  in  model  profiles,  the  effect  of  changing  [03  by  a 
factor  of  two  on  the  wind  calculations  has  been  examined.  In 
Figure  3,  neutral  wind  calculated  with  HSI6  C03  is  compared  to 
that  using  HSIS  C03x2  (which  halves  the  diffusion  velocity).  The 
amplitude  of  the  equatorward  flow  is  damped  (by  up  to  BOm/s  at 
velocity  peaks)  while  that  of  the  poleward  flow  is  slightly 
increased,  but  the  direction  of  flow  (or  the  neutral  wind  pat¬ 
tern)  is  not  significantly  altered.  The  calculated  error  for  the 
experimental  neutral  wind  ranges  from  0  to  30m/s  averaging  around 
16m/s  or  -t-Bm/s.  The  error  bars  are  included  in  the  plots  of 
neutral  winds. 

The  ion  drag  term,  FID,  is  defined  as  the  difference  between 
the  neutral  wind  and  the  ion  velocity  in  the  x  direction  times 
the  collision  frequency,  vni (»Kni (N> ,  where  Kni«BxlOB-8  is  the 
O-t-,0  collision  parameter).  Substituting  for  the  neutral  wind 
from  the  expression  above: 

(4)  FID  «  -CBxl0"®(N)/cosI3 (Vdiff-VzsinI) 

The  local  acceleration  of  the  meridional  neutral  wind,  Uxt,  is 
calculated  by  fitting  a  cubic  to  three  adjacent  measurements. 
Then,  the  net  southward  pressure  gradient  force,  Fx ,  is  calcu¬ 
lated  from  the  momentum  equation  where  the  coriolis  force,  the 
advection  term,  and  the  viscous  drag  are  neglected.  Given  these 
approximations  (discussed  by  Burnside,  19B4) i 


(5) 


Fx  «  Uxt  +  FID 


Nsutral  wind  r*mult«  -from  th»  -four  data  cats  ara  shown  in 
Figura  4.  Soma  ganaral  trands  ara  common  to  all  four  data  sats. 
Maximum  aquatorward  winds  (positiva  Ux)  usually  occur  batwaan 
2300-0100AST.  Staap  accalaration  and  dacalaration  surround  this 
paak.  Winds  drop  off  or  bacoma  pol award  by  0300  followad  by 
southward  accalaration  until  0500.  Tha  daily  minimum  or  most 
polaward  winds  occur  around  noon  (axcapt  for  a  polaward  maximum 
at  6pm  in  January  1985),  otharwisa  oscillating  about  zaro  until 
lata  aftarnoon  when  poleward  flow  occurs  briefly  before  a  rapid 
aquatorward  acceleration  to  the  diurnal  maximum  in  southward 
flow.  This  last  trend  toward  aquatorward  flow  is  interrupted  by 
a  one  hour  poleward  acceleration  at  9pm  which  is  much  more 
noticeable  in  the  winter  data,  causing  100-150m/B  drop  in  equa¬ 
tor  ward  flow  before  the  midnight  maximum  in  aquatorward  flow  is 
reached. 

The  experimental  pressure  gradient  force,  Fx,  is  shown  in 
Figure  5.  In  all  four  data  sets,  the  aquatorward  pressure 
gradient  force  is  strongest  in  the  late  evening  and  the  strongest 
in  the  poleward  direction  around  sunrise  near  equinox  and  in  late 
morning  during  winter  solstice.  (An  aquatorward  (poleward) 
pressure  gradient  force,  Fx ,  implies  a  northward  (southward) 
pressure  gradient.)  A  poleward  gradient  force  is  sometimes 
observed  late  in  the  afternoon,  but  from  noon  to  sunset  the 
gradient  force  is  usually  weak.  The  neutral  wind  exhibits  the 
same  pattern  as  the  pressure  gradient  force  except  during  periods 
of  high  electron  densities  in  the  afternoon  when  the  ion  drag 
force  prevents  the  wind  from  responding  to  the  gradient. 


In  J«nuAry  19B4,  ths  nvutral  Mind  is  •quatorward  at  midnight 
up  to  200in/s  and  bagins  dacalarating  staadily  ravarsing  to  pola- 
ward  flow  at  2am  (saa  Figura  6A) .  On  January  12,  this  poatmid- 
night  abatamant  ia  intarruptad  by  a  auddan  aquatorward  ac- 
calaration  to  lOOm/a  aouthward  flow  at  3am.  Looking  back  at 
Zmax,  it  is  apparant  that  tha  layar  haight  did  not  axhibit  tha 
poatmidnight  drop/riaa  pattarn  found  on  tha  othar  nighta  (Figura 
2B) .  For  tha  othar  January  1984  data,  polaward  flow  occura 
briafly  at  3am  with  maximum  valuaa  of  -'125m/a  bafora  Ux  ravaraaa 
again  and  climba  to  a  aacond  aquatorward  paak  of  150m/a  Juat 
bafora  aunriaa.  In  January  1985,  tha  valocity  at  midnight 
reachaa  325m/B  (400  on  tha  17th)  which  ia  the  diurnal  maximum, 
and  more  day  to  day  variation  ax i eta  in  phaae  and  amplitude 
(Figure  7A) .  The  poatmidnight  poleward  flow  occura  between  3-5am 
and  ia  leaa  than  -25m/a,  the  aecond  aquatorward  paak  ia  lOOm/a 
and  occura  at  5  or  6am,  or  not  at  all.  In  January  1984,  by  8am, 
tha  aquatorward  wind  haa  aubaidad  and  tha  flow  ravaraaa  again 
with  tha  diurnal  maximum  in  poleward  flow  (125m/B)  occurring  at 
noon.  In  January  1985,  the  flow  bacomaa  poleward  between  6-7am 
and  reachaa  a  poleward  maximum  of  •’125m/B  but  returna  to  an 
average  around  zero  by  noon  which  laata  until  5pm  whan  a  auddan 
poleward  acceleration  cauaea  tha  atrongaat  poleward  flow  of  tha 
day,  -175m/e  at  7pm  followed  by  aquatorward  acceleration. 
However,  after  noon  in  1984,  rather  than  a  zero  average  velocity, 
the  poleward  velocity  Juat  dropa  gradually  during  tha  day  and  no 
afternoon  peak  in  poleward  flow  occura  bafora  a  auddan  equator- 
ward  acceleration  which  begina  at  6pm,  rather  than  7pm.  In  both 
January  1984  and  1985,  a  peak  in  aqutorward  wind  occura  at  9pm 
and  ia  alwaya  followed  by  100-150m/a  drop  before  climbing  to  the 
midnight  paak  aquatorward  flow  which  ia  the  diurnal  maximum  in 
1985  (aa  in  the  June  and  September  data  aeta) .  In  1984,  however 
the  atrongaat  daily  winda  occur  with  tha  peak  in  aquatorward  flow 
at  9pm. 

Tha  praaaura  gradient  force  in  January  ia  aouthward  (northward 


prassura  gradiant)  -froni  2ain,  Mhan  tha  wind  ravaraaa  from  pol award 
to  aquatorward  flow,  to  6am.  (January  17,  1985  is  an  axcaption 
in  that  a  strong  poleward  gradient  occurs  at  Sam.)  The  pressure 
gradient  and  ion  drag  forces,  Fx  and  FID,  for  January  1964  and 
1985  are  shown  in  Figures  78, C  and  88, C.  The  resisting  ion  drag 
is  weakly  positive  as  electron  density  is  less  during  this  time 
and  the  wind  responds  to  the  pressure  nradient.  The  effect  of 
electron  density  on  FID  is  apparent  in  comparing  4am  on  January 
17  to  January  19,  1984.  The  pressure  force  is  stronger  on  the 
17th,  but  Ux  is  ISOm/s  on  both  nights.  The  wind  is  not  larger  on 
the  17th  because  the  ion  drag  force  is  a  factor  of  two  larger  on 
the  17th.  (Though  electron  density  at  the  peak  does  not  differ 
significantly  on  the  17th  from  the  19th,  the  perpendicular  ion 
velocity  is  greater  on  the  19th  in  the  north  and  upward  direction 
contributing  to  the  larger  resisting  ion  drag  force.)  The  average 
pressure  gradient  force  in  January  1984  is  zero  by  sunrise  beco¬ 
ming  northward  with  a  -5cm/sec)2(  peak  between  10-1 lam,  an  hour 
before  the  neutral  wind  peaks  in  the  poleward  direction.  Though 
the  pattern  varies  more  rapidly  during  this  time  in  January  1985, 
Fx  is  generally  poleward  with  an  earlier  peak  at  9am.  In  January 

1984,  an  incrase  in  the  ion  drag  between  10-1 lam  reflects  the 
electron  density  which,  as  noted  earlier,  is  highest  in  the 
morning  rather  than  the  afternoon  for  this  data  set.  The  pres¬ 
sure  force  averages  around  zero  in  the  afternoon,  becoming  pole- 
ward  before  reversing  to  equatorward  at  6pm.  In  January  1964, 
the  force  becomes  only  slightly  poleward  before  reversing,  but  in 

1985,  a  very  strong  poleward  Fx  occurs  from  5-6pm  driving  the 
observed  northward  neutral  winds.  As  happens  in  the  morning  in 
January  1984,  the  ion  drag  increases  in  the  afternoon  coincident 
with  maximum  electron  density,  and  when  the  densities  decrease, 
an  increased  FID  reflects  increased  neutral  winds.  In  both  1984 
and  1985,  a  southward  peak  up  to  6cm/e2  occurs  at  8pm  with  the 
force  dropping  off  to  zero  or  even  slightly  negative  by  9pm  and 
climbing  again  to  a  second  southward  peak  force  at  10pm.  The 
wind  exhibits  the  same  pattern  with  roughly  a  one  hour  time  lag. 


A-fter  the  10pn\  peak  in  pressure  gradient,  the  -force  drops  of 
rapidly  preceeding  the  rapid  deceleration  of  Uk . 

Neutral  wind  plots  for  the  September,  19B4,  data  are  shown  in 
Figure  BA.  The  neutral  winds  are  the  strongest  at  midnight  in 
September,  1984,  and  reverse  to  weak  poleward  flow  at  3am,  subsi¬ 
ding  to  zero  and  then  accelerating  poleward  again,  reaching 
lOOm/s  northward  between  6-7am.  As  noted  earlier,  Beptmber  19  is 
a  magnetically  active  day  and  the  winds  undergo  a  sudden  equator — 
ward  acceleration  resulting  in  a  175m/s  southward  velocity 
between  10am  and  noon,  but  the  wind  averages  around  zero  from 
noon  to  sunset,  rather  than  poleward  flow  which  occurs  thoughout 
the  afternoon  in  winter.  Moreover,  on  the  other  September  day, 
after  the  morning  poleward  velocity  subsides  and  reverses  at  9am, 
the  flow  averages  75m/6  equatorward  until  5pm.  The  wind  then 
reverses  and  poleward  flow  <up  to  lOOm/s)  occurs  from  5-9pm. 
This  poleward  reversal  in  the  afternoon  is  also  evident,  though 
less  pronounced,  on  the  magnetically  active  day.  All  September 
data  indicate  equatorward  acceleration  in  the  evening  which 
subsides  at  10pm,  but  increases  again  to  the  strong  equatorward 
flow  (250m/s)  observed  around  midnight  (400m/s  is  observed  on 
September  20) . 

From  Figure  SB  it  is  evident  that  the  experimental  pressure 
gradient  in  September  exhibits  the  same  pattern  as  the  neutral 
winds  from  sunset  to  sunrise,  shifted  an  hour  earlier,  as  in  the 
winter  data.  (The  wind  lags  Fx  by  three  hours  after  midnight, 
September  20  though  the  measured  ion  drag  force,  FID,  is  not 
unusually  large.)  The  maximum  equatorward  pressure  force  occurs 
at  10-1 1pm.  Fx  is  poleward  at  3am,  bam,  and  again  around  bpm.  A 
strong  equatorward  gradient  force  occurs  in  the  afternoon  (3pm) , 
but  does  not  drive  a  strong  acceleration  in  the  neutral  wind 
because  of  high  ion  drag  resulting  from  maximum  electron  density 
at  the  F2  peak. 


Immediately  noticable  in  Figure  9C  is  that  the  June  experi¬ 
mental  neutral  winds  are  almost  always  aquatorward.  Though  there 
is  deceleration  at  sunset,  there  is  not  a  reversal  to  poleward 
wind  before  climbing  to  the  diurnal  maximum  in  equatorward  flow 
(300m/s)  at  midnight.  As  in  the  other  data  sets,  the  wind  falls 
off  rapidly  after  midnight,  but  again,  in  June,  they  do  not 
reverse  at  Sam.  Strong  equatorward  winds  (17Sm/s)  are  observed 
at  6am  followed  by  rapid  deceleration  to  zero  in  the  morning,  but 
not  reversing.  On  one  June  day,  -SOm/s  occurs  around  noon  and  Ux 
averages  slightly  poleward  until  Spm,  but  an  average  50m/s  equa¬ 
torward  occurs  on  the  other  June  afternoon. 

The  experimental  equatorward  pressure  gradient  force  is  often 
a  factor  of  two  larger  in  June  than  in  the  other  data  sets 
(Figure  9B) .  The  diurnal  peak  force  at  10pm  is  10  to  16cm/e  2. 
A  poleward  force  is  dominant  from  midnight  to  sunrise  on  June  26, 
though  the  poet  midnight  period  on  the  other  two  June  data  days 
is  characterized  by  a  weakening  equatorward  force.  Fx  is  weakly 
poleward  from  sunrise  to  noon  on  June  27  concurrent  with  the  only 
poleward  winds  observed  in  June.  Though  an  equatorward  gradient 
persists  past  midnight  on  June  26  and  again  on  June  27,  Ux  begins 
decelerating  at  midnight  due  to  the  high  ion  drag  force 
(FIDslScm/s  2  at  lam  June  27  corresponds  to  unusually  high  peak 
electron  density). 


ThB  nautral  wind  obssrvatians  from  this  study  show  soms  simi- 
laritiss  to  those  from  previous  observations  dicussed  earlier. 
The  magnitudes  are  larger  for  the  current  study.  The  maximum 
neutral  wind  late  in  the  evening  is  usually  around  250m/s  equa- 
torward  while  Harper  found  the  equatorward  velocity  maximum  at 
11pm  to  be  120m/B.  Wind  measurements  from  6t  Bantin-Nancay  (40  N 
lat)  bear  more  resemblance  with  maximum  equatorward  winds  of 
200m/s  in  winter  and  summer  evenings.  Amayenc  found  the  peak 
equatorward  wind  at  equinox  to  be  almost  1/3  less  (120m/s)|  no 
such  seasonal  difference  is  evident  in  the  current  data.  Harper 
found  maximum  poleward  flow  to  occur  at  10am  (October  data)  and 
though  it  occurs  at  this  time  in  the  January  1964  data,  poleward 
flow  in  September  is  strongest  around  sunrise  and  again  sunset. 
Harper  determined  a  flow  reversal  from  equatorward  to  poleward  at 
2am,  a  phenomenon  which  is  observed  at  2-3am  in  this  study. 


Looking  back  at  Figures  2A-D,  it  is  evident  that  the  January  1984 
data  are  unique  in  that  the  maximum  electron  density  of  the  F2 
peak  occurs  Just  before  noon.  The  neutral  wind  and  ion  velocity 
patterns  can  help  explain  this  phenomenon.  The  layer  height  is 
low  from  noon  to  sunset.  More  recombination  is  occur ing  in  the 
afternoon,  reducing  afternoon  peak  densities.  Vertical  ion 
velocity,  Vz  for  January  1984  averages  around  zero  or  weakly 
upward  until  noon  when  downward  motion  begins  and  persists  until 
late  evening.  On  the  other  data  sets,  downward  ion  velocity  is 
not  observed  until  after  6pm.  This  behavior  of  Vz  presumably 
causes  the  observed  lower  Zmax  and  therefore  the  peak  density 
pattern.  The  ion  velocity  perpendicular  to  the  magnetic  field 
lines  can  contribute  to  the  vertical  ion  velocity.  For  this 
January  1984  case,  V-perp  peaks  in  the  north/up  direction  late  in 
the  morning  and  at  noon  begins  to  reverse  to  the  south  and  down¬ 
ward  direction  with  a  -60m/B  peak  velocity  around  3pm  and  an 
average  velocity  in  this  direction  until  after  sunset.  Only  a 
west  electric  field  can  drive  this  component  of  the  plasma 
motion.  The  perpendicular  and  vertical  ion  velocities  for 


24 


January  19,  1984  ara  shown  in  Figura  10  and  ara  rapraaantati va  of 
tha  othar  January  1984  data.  Earliar  obaarvationa  of  daytima 
parioda  at  Aracibo  hava  ahown  a  aouth /downward  V-parp  in  tha 
aftarnoona  implying  a  waat  alactric  fiald  (Behnka  and  Harpar, 
1973).  Only  alactric  fialda  can  driva  iona  parpandicular  to  tha 
magnatic  fiald,  but  tha  nautral  wind  can  contributa  to  tha  vai — 
tical  ion  valocity.  Tha  nc  jtral  wind  in  January  1984  accalarataa 
ataadily  aftar  9ani  to  -125m/a  by  noon.  S/alocitiaa  drop  off 
gradually  but  an  avaraga  pol award  flow  of  -50m/a  laata  through 
moat  of  the  aftsrnoon.  Though  polaward  valocity  raachaa  125m/a 
by  10am  in  January  1985,  the  average  velocity  from  noon  to  sunset 
in  January  1985  is  around  zaro,  and  slightly  aquatorward  in  the 
June  data.  The  poleward  flow  in  the  afternoon  in  January  1984 
pushes  ions  down  the  field  linee  contributing  to  the  downward 
vertical  ion  velocity.  Therefore,  downward  Vz  dominates 
thoughout  the  afternoon  and  recombination  on  the  F2  bottomaide 
occurs  more  rapidly  suppressing  the  afternoon  peak  density. 
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Figure  10.  Vertical  ion  velocity,  Vz,  (positive  up,  cm/s^)  and  perpendicular  ion 
velocity,  V-PERP  (positive  up  and  poleward,  m/s^)  19  Jan  1984. 


COMPARISON  BETWEEN  DATA  AND  MODELS 


Th«  meridional  pressure  gradient  predicted  by  the  MSI8-B3 
model  atmosphere  is  calculated  by  taking  the  north-south  dif¬ 
ference  in  predicted  temperature  and  density.  The  model  force  is 
represented  by  the  solid  line  in  the  plots  of  the  pressure 
gradient  force  for  the  four  data  sets  in  Figures  6Cy7C,6C,  and 
9C.  The  experimental  and  predicted  forces  (proportional  to  minus 
the  pressure  gradient)  are  of  the  same  order  of  magnitude,  in  the 
range  of  ■•'10cm/s2,  though  phases  often  differ.  In  all  four  data 
sets,  both  model  and  observation  indicate  a  peak  in  the  equator- 
ward  pressure  gradient  force  two  or  three  hours  before  midnight, 
though  two  peaks  are  observed  in  the  winter  evening.  The  model 
consistently  predicts  the  most  poleward  gradient  force  (or 
minimum  equatorward  force)  around  6pm  while  it  is  observed  in  the 
late  morning  (except  for  the  afternoon  peak  in  the  poleward  force 
in  January  1985).  MSIS  maintains  a  negative  gradient  force  from 
noon  to  10pm  which  appears  only  briefly  before  sunset  in  the 
data.  The  calculated  Fx  matches  the  model  better  in  phase  during 
the  evening  in  June  and  September  than  in  winter. 


In  January,  the  observed  pressure  gradient  force  has  a  south¬ 
ward  peak  at  45am.  The  MSIS  predicted  pressure  force  peaks  in 
the  southward  direction  at  Tam  (see  Figures  6C,7C) .  Both  obser — 
vation  and  MSIS  indicate  the  force  becoming  more  northward  during 
the  morning.  However,  observed  values  peak  northward  around  10am 
(in  January  19B5  the  pattern  is  varied  but  a  north  pressure  force 
dominates  from  6-lOam)  and  average  around  zero  throughout  the  day 
while  MSIS  indicates  that  the  pressure  force  continues  to 
increase  negatively  to  a  maximum  northward  force  at  6pm.  In  both 
observation  and  MSIS,  the  gradient  force  becomes  southward  in  the 
early  evening.  The  experimental  force  changes  to  southward  at 
7pm  and  quickly  reaches  a  peak  by  6pm  at  which  time  the  model 
force  is  Just  becoming  positive  (southward).  The  observed 
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quickly  ris*s  again  to  a  sacond  southward  paak  at  10pm.  This 
pattarn  is  not  avidant  in  tha  modal  which  continiuas  slowly 
rising  to  a  vary  waak  amplituda  maximum  Just  bafora  midnight. 

In  wintar  tha  pradictad  aquatorward  gradiant  <forca  waakans 
around  sunrisa  to  bacoma  pol award  by  lata  morning,  ramaining 
polaward  until  10  at  night.  Tha  only  polaward  gradiant  forca 
pradictad  in  Juna  and  Saptambar  is  from  3-6  in  tha  aftarnoon.  In 
obsarvation,  Fx  is  polaward  in  Saptambar  bria-fly  from  5-6pm  but, 
in  Juna,  Fx  is  aquatorward  and  incraasing  during  tha  aftarnoon 
and  avening.  Tha  KSI8  force  matches  the  experimental  Fx  in  the 
evening  of  both  June  and  September  in  that  paak  aquatorward  force 
occurs  around  9pm.  liSIS  predicts  a  southward  force  with  minimum 
values  in  the  hours  after  midnight  and  increasing  to  a  maximum 
value  at  Sam  which  in  June,  is  predicted  to  last  until  noon  and, 
in  September,  to  decrease  very  gradually.  After  midnight,  on  two 
of  the  three  June  nights,  the  observed  equatorward  gradient  force 
increases  when  MSIS  predicts  a  minimum  (a  minimum  is  observed  on 
the  other  June  night).  Similarly  in  September,  Fx  has  a  south¬ 
ward  peak  during  tha  modal  prediction  for  a  minimum.  However, a 
smoothed  curve  for  the  September  data  would  have  a  minimum 
(actually  poleward  Fx)  ona  hour  later  than  the  model  minimum  (the 
latter  being  a  weak  equatorward  force).  The  most  marked  dif¬ 
ference  between  the  data  and  the  model  is  that  the  observed  Fx  is 
zero  or  negative  (poleward)  when  MSIS  indicates  a  maximum  posi¬ 
tive  (equatorward)  force  around  sunrise  and  tha  data  begins  to 
increase  positively  as  the  model  becomes  more  negative.  Thus  tha 
phase  of  the  observed  pressure  force  is  opposite  tha  model  from 
sunrise  to  3pm  in  September  and  to  6pm  in  June. 

Dickinson  discusses  NCAR's  Thermospheric  General  Circulation 
Modal,  TGCM,  predictions  of  the  flow  for  100  to  500km  indicating 
circulation  on  constant  pressure  surfaces  at  150  and  300km  (Dic¬ 
kinson,  1984).  Absorption  of  UV  and  EUV  radiation  as  well  as 
auroral  heat  input  are  used  as  energy  sources  with  an  average 


solar  F10.7cm  fluK  of  6B  raprssanting  aquinox  during  solar 
fflinifflum  conditions.  Coincident  geographic  and  geomagnetic  poles 
are  assumed  in  this  version  of  the  model  (denoted  here  as  TGCM- 
1).  At  300km,  neutral  winds  are  shown  to  blow  equatorward  up  to 

BOm/s  from  IBOO  to  0600.  A  slight  northward  component  develops 
after  sunrise  but  drops  off  Just  after  noon.  Acceleration  is 
mostly  in  the  zonal  direction  until  early  evening.  Velocities 
are  more  than  1/3  less  at  150km  blowing  southward  only  from  8pm 
to  midnight.  A  meridional  acceleration  does  not  occur  until  9am 
when  winds  become  northward  up  to  25m/s  dropping  off  by  4pm. 
Where  only  a  slight  poleward  component  to  the  flow  develops  in 
the  morning  at  the  upper  level,  a  dominant  north  acceleration 
persists  throughout  the  day  at  150km.  Conversely,  the  equator- 
ward  acceleration  is  more  persistent  at  300  km.  Figure  11  is 
taken  from  the  published  results  of  this  model. 


Pi9ure  11-  1983  TGCM  prediction  for  equinox  during  solar  mininum.  1200UT  corresponds  to 
08007^T.  Equatorward  flow  is  sore  pronounced  at  300ka;  poleward  at  ISOkn.  (Dickinson,  1984) 


Tim  Killeen  has  analysed  TBCh  including  high  latitude  and  in 
situ  heating  and  assuming  separate  geographic  and  geomagnetic 
poles.  His  analysis  (denoted  here  as  T6CM-2)  determines  among 
other  parameters,  predicted  neutral  winds  specifically  for  Are- 
cibo.  Another  generation  of  TBCH  by  Casandra  Fesen  (referred  to 
here  as  T6CI1-3)  produces  vertical  cross  sections  of  velocity 
fields  for  Arecibo  assuming  coincident  geographic  and  geomagnetic 
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pol*s,  •xcluding  high  latitude  heating*  and  excluding  or  inclu¬ 
ding  the  e-f<fects  o-f  tidal  propogation  4roin  below  the  thermosphere 
(the  most  recent  model  development).  Model  results  are  shown  in 
Figure  12A  -for  solar  maximum  equinox.  The  experimental  winds  are 
in  geomagnetic  coordinates.  The  geomagnetic  coordinates  are 
plotted  to  emphasize  differences  in  the  model  run  with  separate 
geographic  and  geomagnetic  poles  and  that  assuming  coincident 
poles  (for  which  geographic  coordinates  are  plotted).  Difference 
in  these  two  plots  also  arises  from  different  type  graphs  used  to 
draw  each  and  to  a  less  extent  from  excluding  high  latitude 
heating  from  the  TGCM-3  analysis.  At  SOOkm*  in  geomagnetic 
coordinates,  poleward  flow  is  indicated  during  the  day  from  7am 
to  2pm  with  a  45m/s  peak  at  9am.  The  flow  is  southward 
throughout  the  rest  of  the  24  hour  period  reaching  a  maximum  of 
45m/s  at  3am.  The  geographic  poleward  flow  begins  an  hour 
earlier  and  lasts  four  hours  longer  in  the  afternoon  and  the 
magnitude  of  the  peak  equatorward  velocity  is  reduced  by  a  half. 
One  might  keep  in  mind  these  differences  in  comparing  observation 
to  the  remaining  model  plots  which  are  all  in  geographic  cordi- 
nates.  The  tidal  effects  contribute  to  the  duration  of  the 
poleward  meridional  flow,  shifting  the  maximum  poleward  flow  to 
the  afternoon  rather  then  in  the  late  morning  but  with  little 
change  in  magnitude.  The  magnitude  of  the  equatorward  velocity 
nearly  doubles  from  that  without  tide  effects,  but  the  duration 
of  the  equatorward  flow  is  reduced,  lasting  only  from  7pm  to  3am. 


Various  T6CM  results  for  solar  minimum  equinox  are  shown  in 
Figure  12B.  The  TGCM-1  prediction  is  derived  from  the  global  map 
in  Figure  11  by  taking  longitude  to  be  equivalent  to  a  time  scale 
and  measuring  the  meridional  wind  component  over  Arecibo. 
Another  plot  in  Figure  126  was  drawn  from  a  vertical  cross 
section  produced  from  TGCM-3  without  tides.  The  trend  in  the 
neutral  wind  from  these  model  results  is  for  poleward  flow  begin¬ 
ning  at  4am,  reaching  a  maximum  in  late  morning  ('^45m/s)  and 
decreasing  in  the  afternoon  reversing  to  equatorward  around  5pm. 
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Th*  predicted  peak  equatorward  velocity  is  at  midnight,  and 
though  the  recent  model  only  indicates  peak  epeeds  of  25m/s,  in 
geomagnetic  coordinates,  the  predicted  maximum  would  be  greater. 
Also  included  in  Figure  12B  is  a  plot  drawn  from  the  vertical 
cross  section  at  Arecibo  from  T6CM-3  wiht  the  effects  of  tides 
propagating  from  below  included.  The  poleward  flow  begins 
earlier  (2am)  and  peaks  at  4am  and  5pm  with  magnitude  comparable 
to  that  for  the  peak  poleward  flow  predicted  to  occur  in  the  late 
morning  by  modeling  without  tides.  The  equatorward  velocities 
are  doubled  with  the  peak  predicted  to  occur  an  hour  before 
midnight. 

The  experimental  neutral  winds  for  the  four  September,  19B4 
days  characterize  solar  minimum  equinox  for  a  comparison  between 
theory  and  observation  (Figure  136).  (The  model  uses  F10.7bB0 
and  the  September  observations  are  for  FIO. 7^100.)  The  phase  of 
the  September  winds  matches  the  theoretical  plot  which  includes 
tidal  effects  (solid  line  in  Figure  12S) |  however,  observed 
amplitudes  are  a  factor  of  three  larger.  The  maximum  predicted 
equatorward  winds  are  6Sm/s  while  peak  observed  winds  are  usually 
250m/s  (400m/s  has  been  observed) .  The  strongest  equatorward 
flow  is  at  midnight  in  both  theory  and  observation,  though  the 
theory  calls  for  peak  poleward  flow  at  3am  and  though  poleward 
flow  is  observed  at  that  time,  stronger  poleward  flow  is  observed 
at  7am. 

The  T6CM-2  prediction  (with  the  high  latitude  heating  and 
separate  poles)  has  also  been  analysed  for  Arecibo  for  winter 
solstice  solar  maximum  conditions  (Figure  120.  The  model  indi¬ 
cates  poleward  flow  all  day  and  night  except  for  two  hours  around 
3am  when  the  meridional  wind  drops  to  zero  or  lOm/s  southward. 

The  observations  of  neutral  wind  in  January  1984,  Just  fol¬ 
lowing  winter  solstice  and  two  and  a  half  years  after  solar 
maximum,  show  little  similarity  in  phase  to  theory,  and  model 
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winds  sre  always  polaward,  while  except  for  a  brief  post-midnight 
reversal,  the  observed  winds  are  equatorward  from  sunset  to 
sunrise  (Figure  13A) .  The  peak  in  equatorward  winds  observed  at 
6am  could  correspond  to  the  theoretical  minimum  poleward  flow 
indicated  two  hours  earlier.  Similarly,  a  maximum  in  poleward 
flow  is  predicted  at  Bam,  but  observed  at  noon.  A  maximum  in 
poleward  wind  is  evident  at  6pm  in  both  theory  and  observation. 
Though  the  model  poleward  wind  decelerates  in  the  evening,  theory 
does  not  indicate  the  maximum  equatorward  wind  observed  at  6pm 
and  again  at  midnight. 
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Figure  12A.  Dotted:  TGai-2  (geon.  coorda.) 

Hashed:  TGCM-3,  without  tldee 
solid:  TCCM-3,  with  tldee 


Fur  rcfprfi:co  to  Figures  12A-C  (also  see  text): 

TOttl-l:  (Dickinson)  Model  lnclU(1>'S  higli  latitui)<’ 
and  In  eltu  heating,  coincident  geog/qeom 
|K>les.  Plot  drawn  from  global  map  in 
Figure  11. 

TGCM-2:  (Killeen)  Model  includes  hlgli  latitude 
and  in  situ  heating,  separate  geog/geom 
poles.  Plots  drawn  from  graphs  of  Ux  vs 
time  at  Arecibo. 

TGC)4-3:  (Fesen)  Model  Includes  in  situ  heating, 
coincident  geog/qeom  poles,  with/without 
tides.  Plots  drawn  from  vertical  cross 
sections  of  contours  of  Ux  in  time. 
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Figure  12B.  Dotted:  TGCM-2 

Dashed;  TGCM-3,  without  tides 
Solid;  TGCM-3,  with  tides 
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Figure  13A.  January  1984  observations  to  diaracterize  winter  solstice  two  years 
before  solar  sUnimum.  The  model  results  are  for  winter  solstice  solar  maximum  and 
do  not  include  effects  of  tides  propoqating  from  below  (Pig  12C)  and  show  little 
similarity  to  these  observations. 
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Figure  13B.  Septenber  1984  observations  ciiaracteriec  equinox  solar  minimum  and 
can  be  compared  to  the  solid  line  in  Fig  12b.  The  pliase  in  theory  is  similar  to 
observation,  but  the  observed  amt'litude  is  a  factor  of  three  larger  and  there  is 
an  equatorward  shift  of  the  serolinc. 
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THEORY 


Th«  tcndsncy  toward  aquatorward  winds  at  Aracibo  in  ths  present 
data  might  be  related  to  a  Hadley  cell  circulation  driven  by  high 
latitude  heating  in  which  a  return  poleward  flow  occurs  at  lower 
levels.  The  previously  mentioned  wind  contours  for  Arecibo 
during  solar  maximum  showed  predominantly  equatorward  flow  above 
and  poleward  flow  below  300km  (Burnside  et  al . «  1983).  Global 
maps  from  TGCH  for  ths  constant  pressure  surfaces  at  150  and 
300km  also  show  more  pronounced  northward  flow  at  the  lower  level 
with  more  southward  flow  aloft  (Dickinson,  1983).  Hayr  discusses 
a  Hadley  cell  circulation  with  the  upper  branch  having  poleward 
flow  (Mayr  et  al . ,  1984).  He  suggests  that  with  solar  dif¬ 
ferential  heating  at  low  latitudes,  temperature,  and  pressure 
decrease  toward  the  poles  causing  a  Hadley  cell.  The  redistri¬ 
bution  resulting  from  wind  induced  diffusion  depletes  the  atomic 
oxygen  concentration  at  ths  low  latitudes  thus  possibly  reversing 
the  pressure  gradient  (cf.Mayr  et  al . ,  1984). 

Early  determinations  of  thermospheric  wind  speed  from  analysis 
of  changes  in  satellite  orbit  inclinations  indicated  a  mean 
eastward  zonal  wind  of  about  130m/s  at  300km  for  30  latitude 
(King-Hsle  et  al .  ,  1970).  This  would  result  in  a  component  of 
18m/s  in  the  magnetic  south  direction  over  Arecibo,  but  is  not 
enough  to  explain  differences  in  TGCH  predictions  and  the  current 
observations.  In  addition,  theoretical  models  produce  much 
smaller  values  of  mean  zonal  winds  supported  by  mass  spectrometer 
measurements  from  Dynamics  Explorer  II  (Hayr  et  al . ,  1984). 

The  observed  temperatures  agree  well  with  HSIS  in  both  phase 
and  amplitude.  However,  although  the  observed  pressure  force  is 
of  the  same  magnitude  as  that  predicted  by  HSIS,  the  phases  are 
significantly  different.  Rishbeth  (1977)  determined  that  a  more 
complex  pressure  distribution  than  ths  Jacchia  model  then  being 
used  would  be  necessary  to  explain  Arecibo  nightime  winds. 


Parhaps  (ISIS  can  not  raproduc*  small  seals  dynamics  influsncing 
the  equatorial  tropical  wind  -field.  Rishbeth  describes  a  tro¬ 
pical  high,  a  diurnal  bulge  in  the  pressure  -field  -from  which 
diverging  wind  would  give  a  poleward  -flow  at  Arecibo.  He  sug¬ 
gests  that  a  displacement  of  this  high  would  result  in  transequa- 
torial  flow  and«  therefore,  for  a  high  situated  north  of  Arecibo, 
an  equator ward  flow  component. 

The  T6CM  underestimates  the  amplitude  of  the  neutral  wind 
observed  in  this  study.  The  inclusion  of  high  latitude  heating 
and  separate  geomagnetic  and  geographic  poles  in  the  model  used 
for  comparison  would  increase  the  predicted  amplitude  some,  but 
not  nearly  by  a  factor  of  three.  The  neutral  winds  in  this  study 
are  deduced  from  the  measured  parallel  ion  velocity  and  the 
calculated  diffusion  velocity.  The  magnitude  V-par  is  consistent 
with  previous  results  <<100m/s)}  the  diffusion  velocity  is 
generating  the  large  neutral  winds  observed.  The  sensitivity  of 
the  neutral  wind  to  atomic  oxygen  density  has  been  demonstrated. 
If  the  densities  predicted  by  MSIS  are  too  small,  through  the 
diffusion  coeficient,  V-Diff  and  therefore  Uk  will  be  too  large, 
though  again,  not  by  a  factor  of  three.  On  the  other  hand,  the 
TGCn  model  could  be  overestimating  electron  density  and  therefore 
overestimating  ion  drag  resulting  in  smaller  winds.  In  addition 
to  bringing  the  phase  in  close  agreement  with  observation,  the 
inclusion  of  tides  propagating  from  below  the  thermosphere 
increases  the  neutral  wind  amplitude.  The  forcing  term  for  the 
semidiurnal  tidal  modes  used  is  tuned  to  observations  at  Arecibo 
which,  as  noted  previously,  are  much  less  than  those  observed  in 
this  study. 
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